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Abstract 

We study the time dependent polariton condensation as well as the parametric scattering process of polaritons in a semiconductor 
microcavity. Based upon a new stochastic scheme the dynamics for both cases is fully analyzed. We show how the evolution of 
the system is described by a set of stochastic differential Schrodinger equations which in average reproduces the exact dynamics. 
Furthermore, we underline the role that Coulomb correlations plays in the polariton dynamics. Threshold behaviors are well 
captured by the present approach. The results are in complete agreement with recent experimental observations. 
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1. Introduction 

Microcavity polaritons are quasiparticles created in the 
strong coupling regime between photons in a microcavity 
and the excitonic resonance of a semiconductor quantum 
well [1]. Polaritons satisfy Bose statistics (at least in the 
low density regime) and are excellent candidates to study 
quantum effects at the macroscopic level, e.g. Bosc-Einstein 
Condensation (BEC) and superfluidity. As a practical ad- 
vantage that polariton systems exhibit over other quasi- 
bosonic particles in condensed matter systems is their easy 
experimental optical control [2,3]. Particulary, polaritons 
can be optically injected in the microcavity by external 
light impinging on the semiconductor nanostructure, and 
their properties can be inferred from the observation of the 
emitted light. The possibility of obtaining a quantum fluid 
in a solid state system, with easy control and integration 
potentialities, could also open a new promising way to the 
implementation of quantum information technologies. 

The polariton condensation has being actively pursued, 
since no BEC has being observed for a solid state system. 
Recently, a possible quantum phase transition of polaritons 
was suggested [4] , and the spontaneous formation of a co- 
herent polariton state was reported in [5]. This last result 
shows a remarkable manifestation of the bosonic behavior 
of this short living particles, but it has been interpreted 
as a polariton laser phase transition. From the theoretical 
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side, most existing works [4] pursue a strict analogy with 
the laser theory, where the quantum fluctuations responsi- 
ble for the phase transition appear directly from the parti- 
cle number conservation. 

On the other hand, it is well known that due to their 
excitonic component, polaritons are subjected to Coulomb 
interactions. This nonlinearity is expected to trigger the 
quantum phase transition. Moreover, the unexpected ob- 
servation of a thermal-like intensity far above the threshold 
[5] does not match with the laser picture, leaving the unan- 
swered question of the role of interactions, usually over- 
looked in the theory. A more recent attempt to address 
this question has been made in [6] , where the laser picture 
was left behind and polariton-polariton interactions were 
taken into account. Nonetheless, this last study is within 
the scope of a mean- field approximation (MFA), and thus 
fails to give a proper description of the polariton dynamics 
above the parametric threshold where high order Coulomb 
correlations are more important than in the low density 
regime. 

The present paper is devoted to gain more physical in- 
sight into the polariton dynamics by analyzing two differ- 
ent dynamical processes: (i) the time-dependent condensa- 
tion and coherence build-up of polaritons relaxing into the 
lowest branch states and (ii) the time-dependent paramet- 
ric scattering. Both processes will be considered within the 
same theoretical approach based on the stochastic unrav- 
elling of the dynamics of an interacting many-body system 
embedded in a noisy environment. Partial realizations of 
such a computational scheme have been previously consid- 



ered. On the one hand, closed quantum interacting systems 
have been studied by stochastic unravellings in Ref.[7] for 
bosons and in Rcf.[8] for fermions. On the other hand, sim- 
ple open quantum systems have been studied by following 
the stochastic density matrix for the whole system plus 
environment in Rcfs.[10,ll]. It is one of the main aims of 
the present work to extend those partial developments to 
a complex many-body open quantum system such as the 
one offered by microcavity polaritons. In order to provide 
quantitative predictions, previous stochastic schemes are 
extended to follow the evolution of a polariton system cou- 
pled to a noisy environment where intcrparticle interactions 
and environment effects can be treated on the same foot- 
ings. This will allow us to transform the evolution of an open 
interacting polariton system, into the stochastic evolution 
of a single polariton [7] separately from the stochastic evo- 
lution of the environment [10,11]. As compared with simi- 
lar approaches the advantages of this stochastic scheme are 
manifold: first of all, the complete dynamics, both transient 
and steady-state regimes, can be obtained with full fluc- 
tuations due to interactions properly taken into account. 
Second, a many-polariton problem can be mapped onto a 
single polariton system driven by stochastic terms. 

The paper is organized as follows: in Section 2 a gen- 
eral description of the stochastic scheme is presented. In 
Section 3 the polariton condensation dynamics with both 
Coulomb interactions and phonon scattering is considered. 
We pay particular attention to the quantum properties of 
the condensate fraction and the second-order degree of co- 
herence. The parametric scattering of pump polaritons in 
the lower polariton branch and the dynamics of formation 
of the eight-shaped final allowed states is studied in Section 
4, where comparison to experimental results [2] reveals the 
importance of including the nonlincarities associated to po- 
lariton interactions, specially above the threshold regime. 
Finally, in Section 5 we sketch our concluding remarks. 

2. Stochastic scheme for interacting bosons 

We start with a brief review of the main assumptions for 
a stochastic scheme aiming to determine the evolution of a 
closed many-boson system [7]. The initial full information 
of the system is encoded in the state vector — \N : (j>) 
where N bosons arc accommodated in the same single par- 
ticle state \<p). Thus, one needs only to find the appropriate 
evolution of \4>) in order to obtain the N-bosons system dy- 
namics. For non interacting particles this is of course triv- 
ial, but in the interacting case the action of the Hamilto- 
nian on the N-particlc state vector produces terms orthog- 
onal to \4>). This cannot be taken into account by standard 
MFA, and hence the necessity to introduce a noise term 
with specific properties. The first non-trivial application of 
this method was performed in Ref.[12]. 

For microcavity polariton systems at low temperatures, 
an always present source of nonlincarity is the coupling 
with the acoustic phonon reservoir. At first glance, the ne- 



cessity to include this reservoir, should leave the stochastic 
scheme just mentioned out of the picture for polariton dy- 
namics. However, as it has recently been stated for simple 
open systems [10,11], the quantum system and its environ- 
ment can be described by a pair of decoupled Stochastic 
Schrodingcr Differential Equations (SSDE). The success of 
this fictitious evolution is based on the adequate introduc- 
tion of stochastic terms. This will guarantee that when av- 
eraged over the stochastic trajectories of both subsystems, 
the exact whole dynamics is recovered. The central quan- 
tum system, however decoupled from its surrounding in 
every individual trajectory, will feel the effect of the envi- 
ronment through a driven noise term, allowing in this way 
to obtain the exact reduced density matrix at any time. In 
the following of this Section we will give a brief guideline 
on how the two just mentioned stochastic schemes can be 
extended to consider polariton dynamics. 

In general, the full Hamiltonian can be written as 

H T = H S +H E + Hi (1) 

where H$ describes the central quantum subsystem (in- 
cluding intcrparticle interactions), He is the Hamiltonian 
associated to the environment and Hi is the coupling 
Hamiltonian which we assume can be expressed as Hi = 
J2 a ® B a , where A and B are operators acting on the 
system and environment, respectively. In order to obtain 
the formally exact density matrix as the stochastic average 
over multiple realizations, each sample has to be taken as 
the dyadics p = |$i)($2|- Moreover, one of the main ad- 
vantages of the present approach is that the different kets 
and bras in the dyadics can be taken as the separable forms 
\<&v) = IVv) ® \Xv) {v = 1,2) of system and environment 
state vectors, respectively. Of course, the stochastic aver- 
age reproduces the true entanglement between quantum 
system and environment, as it should be. The evolution of 
the whole system under the Hamiltonian in Eq.(l) is then 
given by the set of SSDEs [10] 

d\Mt)) = -idtH s \i/> v ) -idtJ2dW a „(t)A a \Mt)) (2) 

a 

d\Xu{t)) = -idtH E \ Xv ) + y £ J dW av (t)B a \ x „{t)) (3) 

a 

with noise terms satisfying E[dW a i v < dW au (t)] = 8 a ',a&u',u 
and E[dW ay ] = (E[...} means the average over stochastic 
realizations) . By averaging the dyadics built on the states 
evolving as in Eq.(2) and Eq.(3), the exact Liouville equa- 
tion for the whole density operator is found. This has been 
extensively discussed in Ref.[10] within a quantum jump 
approach and alternatively in Ref. [11] using a diffusion 
method. 

Assuming the Hamiltonian H$ in Eq.(l) conserves the 
number of particles, one can consider the dynamics of a 
closed quantum system due to the separability from the en- 
vironment. Hence, the new idea is to describe the evolution 
of the N interacting particles, now in an open situation, as 
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occupying the single particle state vector \<f>). In order to 
find the stochastic evolution within this prescription, the 
new form of the state vector \ip) = \N : <j>) has to be adapted 
to Eq.(2), verifying that the SSDEs once again yield to the 
formally exact dynamics. This procedure implies that 

d\<j>„) = -idtG v \^ v ) + dM v \(j) v ) - iJ2a a \^ v )dW av (t) (4) 

a 

where G v contains inter-particle interaction effects treated 
at the level of MFA while dM v is a noise term responsible for 
capturing the particle-particle interaction effects beyond 
MFA. Finally, the third term in Eq.(4) carries the net effect 
of the environment coupled to the system. The notation a a 
is used to represent a system's coupling operator projected 
into the single particle space. 

It is important to note that within this stochastic ap- 
proach, where system and environment are always separa- 
ble for each individual realization, the calculation of any 
system's observable, after tracing over the environment, is 
given by 

{6(t)) = E[(Mt)\d(t)\Mt))(x2\xi)}- (5) 

With this in mind, a complex problem where bosonic par- 
ticles interact with each other and with their surroundings, 
has been reduced to a set of SSDEs for a single particle, 
Eq.(4), decoupled from its environment, Eq.(3). The price 
to be paid for this formidable reduction of complexity is 
that a large number of stochastic realizations has to be per- 
formed in order to recover the exact solution. Obviously, the 
results obtained by this procedure are beyond MFA while 
additionally particle-particle interactions as well as envi- 
ronment effects (phonons) can be treated on equal grounds. 

3. Spontaneous coherence buildup of microcavity 
polaritons 

In a typical photoluminesccncc experiment under non- 
resonant excitation, the strong energy dispersion of micro- 
cavity polaritons leads to a bottleneck relaxation dynamics 
[14] , where polaritons go to the flat exciton-like region of the 
energy dispersion curve. From this point on, polaritons re- 
lax into the bottom states of the lower branch. Using these 
considerations a simplified model of polariton dynamics is 
an effective two- level system [3], where the ground state 
corresponds to k = which is non-degenerate {Eq), and 
the excited or bottleneck set of states is macroscopically 
degenerate {E{) (see Figure 1). 

In the lower polartion branch the Hamiltonian describing 
both Coulomb and phonon scattering is 

H = Ukp{pk + u i b \ h i + H c+ H ph (6) 

fe q 

H c = \Y. kk'qv^, P l +q p{,_ qPk 'P k (7) 
H ph = Y,9 ( k i{b\ P W+H.C.) (8) 

kk' q 
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Fig. 1. Lower polariton energy branch versus total momentum in a 
linear-log scale. The excited state is represented by the effective bot- 
tleneck state, which is at the edge of the flat exciton-like dispersion 
region. The ground state is the zero momentum state in the bottom 
of the dispersion curve. 

where p k is the Bose annihilation operator for a polariton 
with momentum k. The excitonic component of the polari- 
tons is responsible for Coulomb interactions, which is cap- 
tured in H c with the coupling strength v^ k \ . The semicon- 
ductor environment is modelled by acoustical phonons dis- 
tributed according to the lattice temperature. The corre- 
sponding phonon annihilation operator is b q with the usual 
Bose commutation relations [b q , b ql ] = 5 qq / . The polariton- 
phonon scattering is responsible for decoherence and is rep- 
resented in Eq.(6) as H v n with g^, as the coupling strength. 
Note that this polariton Hamiltonian is number conserv- 
ing. The dynamical condensation of microcavity polaritons 
will be studied using the stochastic formalism exposed in 
the latter Section. The relevant material parameters [2] are 
typical for a GaAs-based microcavity: heavy-hole exciton 
Rabi splitting HQr = l.&2meV, the detunning with the 
confined photon mode —2AmeV and the effective index of 
the cavity mode 3.5. The quantization area A = 100/^to 2 
determines the Coulomb interaction [13] and the phonon 
coupling terms [14] . 

We consider N interacting polaritons in the presence of 
a semiconductor environment at zero temperature. Since, 
the effective two level model does not consider the relax- 
ation from high energy states into the effective bottleneck 
state, the initial state consists of the total number of po- 
laritons randomly distributed in the excited energy level 
(importantly, no seed is needed for level Eq). Indeed, a 
zero temperature environment will force the polaritons into 
their ground state following a non-trivial evolution due to 
polariton interactions. Our main interest is to give a full 
description of the ground state statistics. We first calculate 
the ground state probability distribution Po(n^t). This is 
initially a delta function centered at n = 0, since all polari- 
tons are initially in the excited state. For comparison pur- 
poses, an overall factor, 7, is introduced in the Coulomb in- 
teraction term in order to account for the effective strength 
of two body Coulomb interactions: 7 = corresponds to 
no Coulomb interaction. As is clearly shown in Figure 2, 
for a fixed total number of polaritons N = 1600, the po- 
lariton system decays monotonically into its ground state 
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Fig. 2. Ground state probability distribution Po(n, t) for different 
polariton interaction strengths: (a) 7 = 0, (b) 7 = 0.4, (c) 7 = 0.8 
and (d) 7 = 1, averaged over 5 X 10 5 realizations. The total number 
of polaritons is N = 1600 and lattice temperature T=0. 
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Fig. 3. Projections of the ground state probability distribution 
Po (n,t) at four various times: 100 ps (green), 200 ps (red), 300 ps 
(blue) and 400 ps (purple). Results are shown for N=1600, 5 X 10 5 
realizations and different polariton interaction strengths: (a) 7 = 0, 
(b) 7 = 0.4, (c) 7 = 0.8 and (d) 7 = 1. 

when interactions are not included. As the interactions are 
turned on, the system's dynamics changes completely. Po- 
laritons start with a slow relaxation dynamics, but after 
a finite time the ground state becomes rapidly populated 
with a macroscopic fraction of N. 

The Gaussian-like distribution for vanishing inter- 
particle interactions is evidently different from the Poisson- 
likc distribution when the interactions are fully included, 
as displayed in Fig. 3. This latter feature is a clear signa- 
ture of an interaction driven coherence buildup process for 
which an off-diagonal-long-range order (ODLRO) can de- 
velop. It is well known [15] that ODLRO does not implies 
a finite value of (po)i m agreement with our results. 

The so-called coherence degree parameter is defined as 
i] = 2 - 5 (2) (0), with 5 (2) (0) = {p^pD- This parame- 
ter takes values between and 1 for states ranging from 
chaotic/ thermal to coherent. Using this parameter the full 
statistical properties of the ground state arc unravelled. 
Starting with an incoherent state (N polaritons randomly 



Fig. 4. Time evolution of the coherence degree parameter r): (a) 
N = 1600 and different interaction strengths 7 = (red), 7 = 0.4 
(green), 7 = 0.8 (blue) and 7 = 1 (purple), (b) Different number of 
polaritons (7 = 1) N = 1400 (red), 1600 (green) and 1800 (blue). 
Lattice temperature T=0 and 5 X 10 5 stochastic trajectories. Dis- 
persion bars are shown for the intermediate case of N = 1600. 

distributed in the excited level) , the coherence rises as the 
ground state becomes populated. This is clearly shown in 
Fig. 4a. Again, the most interesting observation about this 
evolution is the difference between the smooth growing of 
coherence in the non-interacting regime as compared with 
the sudden rise of coherence for the interacting case, com- 
plementing the results depicted in Fig. 3. A threshold be- 
havior is clearly observed in Fig. 4b. Once the number of po- 
laritons, or cquivalently the density, is increased beyond a 
certain value the coherence rises suddenly A possible phase 
transition triggered by Coulomb interactions is thus rein- 
forced. Moreover, since no high order correlations are in 
principle discarded by the present method, as it is usually 
the case for condensed matter systems [16,17] (excitons) 
and [6] (polaritons), we can quantify their influence on the 
polariton condensation dynamics. 

Finally, we would like to remark some of the most im- 
portant aspects of these results. First of all, the threshold 
behavior that governs the dynamics was already discussed 
in Ref. [4] , and reported results are in good agreement with 
ours. This gives us confidence that the stochastic scheme we 
used is indeed a very practical tool for the study of complex 
systems like those found in semiconductor microcavities. In 
particular, effects due to the semiconductor environment 
are responsible for the relaxation of the system, leading to 
a continuous accumulation of polaritons in the same final 
state at low temperatures. Nevertheless, correlations gen- 
erated by interactions play a crucial role in the polariton 
condensation dynamics. While a coherence degree parame- 
ter close to 1 can only be interpreted as a criteria for polari- 
ton lasing [4] , instead of BEC in a quasi-two-dimensional 
nanostructure, signs of a true phase transition are present 
in our results as a consequence not only of the bosonic be- 
havior but from quantum fluctuations provided by polari- 
ton interactions. Since the appearance of the phase transi- 
tion is governed by Coulomb interactions, our results point 
toward a BEC type behavior rather than a polariton laser. 



4. Polariton parametric scattering 

In recent experiments [2,3] polaritons are optically ex- 
cited at a well defined in-plane momentum k, allowing a 
direct control of their dynamics. The strong Coulomb in- 
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Fig. 5. Contour plot of the energy balance equation 
| + Ek i — 2i?fe p | . The blue line represents the set of final states 
for two pump polaritons at k p = {0.66, 0}fim~ 1 spontaneously scat- 
tered conserving both energy and momentum. 

teraction results in a spontaneous parametric scattering of 
polariton pairs (signal and idler) according to momentum 
and energy conservation criteria. For resonant optical ex- 
citation, the scattering leads into a final state 8-shaped 
distribution in k space, as first suggested by Ciuti et al 
[13] and experimentally confirmed in Ref.[2]. We use the 
present stochastic scheme to study the dynamics of initially 
pumped polaritons in a particular momentum state k in 
the lower polariton branch. Once again, the advantages of 
the present approach are manifold. Similar to the conden- 
sation of polaritons studied in the previous Section, the the- 
oretical discussion of parametric scattering has been based 
upon MFA, and thus the results are only valid below the 
threshold of parametric luminescence. Thus, the study of 
the macroscopic behavior of interacting polaritons in a re- 
laxing environment has been poorly investigated [9] . More- 
over, the dynamical behavior of condensed matter quasi- 
particles in an out of equilibrium regime can be systemati- 
cally studied within our approach. In this Section we aban- 
don the effective two-level approximation and resort to an 
extended multi-level description. 

For the numerical simulations that follows, the total 
Hamiltonian is still given by Eq.6, where the system's 
parameters are the same as in Section 3. The possible 
final states for the scattering of two polaritons pumped 
at k = {fc p ,0}, into a signal and idler pair, are displayed 
in Fig. 5, where the long-time limit, in which both energy 
and lineal momentum have to be conserved, is depicted 
as the continuous line. This information is crucial for un- 
derstanding the many-polariton case that we discuss from 
now on. The initial state is taken as N polaritons pumped 
at k p = {0.66, 0}/im _1 . The system is let to evolve and 
population in k space is monitored at different times. In 
Fig. 6 contour plots of the polariton population distribution 
are shown, for N=l 000 particles initially pumped. Clearly, 
a strong asymmetry in signal and idler emission patterns 
is seen. In particular, although there is a small population 
in the idler modes, it is not comparable to the signal one. 

Above the parametric threshold a remarkable different 
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Fig. 6. Population evolution in momentum space of N=1000 polari- 
tons initially pumped at k = {0.66, 0}/im _1 . The plots correspond 
to different times: (a) 20 ps, (b) 40 ps, (c) 60 ps (d) 80 ps, and 10 4 
simulations. 



behavior is displayed. As the density of polaritons increases 
a more significant fraction of the population is concentrated 
in lowest momentum states, and thus an important inten- 
sity is registered in the idler state opposed to the signal 
one, as can be seen in Fig. 7. Similar results were reported 
in Rcf. [9] . Furthermore, the finite intensity that rises in the 
idler state is asyncronous with the evolution of the signal 
population, in close agreement with experimental results in 
Rcf. [2]. The asymmetry in the shape and time of the for- 
mation dynamics of the final state distribution in k space is 
a direct result of the different multiple scattering processes 
affecting signal and idler polaritons. As is well known, for 
higher momentum states (idler modes), polaritons have a 
big excitonic component, which implies a stronger coupling 
with the lattice phonons, and as a consequence a fast de- 
cay in idler population. Even though in the time range of 
the reported dynamics the dominant scattering processes 
correspond to Coulomb interactions, some evidence of de- 
coherence in the idler states, due to the semiconductor en- 
vironment (phonons), is visible in Fig. 7c and 7d. 

In contrast to MFA, the threshold behavior is well cap- 
tured by our stochastic method. As the polariton density 
increases, a dramatic change is observed in the signal pop- 
ulation. While below the threshold the signal population is 
weak, a sudden increase close to the threshold is observed. 
In Fig. 8, the signal population is displayed for different 
times as a function of the number of initially pumped po- 
laritons. Close to N — 2000, a clear change of slope is reg- 
istered for any time. The behavior of the idler mode popu- 
lations (not shown) is analogous. 
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Fig. 7. Population evolution in the momentum space of N=5000 
polaritons initially pumped at k = {0.66, 0}/im _1 . The plots corre- 
spond to different times (a) 20 ps, (b) 40 ps, (c) 60 ps, (d) 80 ps, 
and 10 4 simulations. 



many body problems in condensed matter scenarios. Its 
validity was tested and compared positively with previous 
theoretical and experimental results. 

In the dynamical condensation of polaritons, we found 
that while the semiconductor environment induces a slow 
relaxation process, polariton-polariton interactions trigger 
the sudden development of a macroscopic coherent state. 
This differs clearly from the polariton laser picture, and 
other recent proposals which claim that this behavior is 
just a consequence of the particle number conservation [4]. 
While a more detailed study on this phenomenon is needed, 
our results favor the BEC type phase transition of micro- 
cavity polaritons. 

On the other hand, the dynamical process of paramet- 
ric emission was studied. There exists an asymmetry in 
shape and time dependence of the far field emission of sig- 
nal and idler modes governed by this kind of scattering 
processes. The threshold behavior is also well captured by 
our stochastic method. This corresponds to the sudden rise 
in signal and idler populations. Coherent properties of this 
parametric emission process may also be studied by the 
present method. Additionally, a straightforward possibility 
is to extend the present stochastic method to finite lattice 
temperature situations. 

This work was supported by COLCIENCIAS-1204-05- 
13614 (Colombia) and Facultad de Cicncias-UniAndcs. 
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Fig. 8. Population fraction of the signal mode k s = {— 0.5, 0}/im _1 
as a function of the initially pumped polaritons N at different times: 
t = 5 ps (gray), t = 10 ps (green), t = 15 ps (blue) and i = 20 ps 
(red). f0 4 realizations for each initial number of total polaritons. 

5. Conclusions 

In the present paper we have presented a stochastic study 
on the dynamical properties of microcavity polaritons. In 
particular, we reviewed the dynamical condensation of po- 
laritons emphasizing on the role of Coulomb interactions. 
Additionally, the threshold behavior of the spontaneous 
parametric downconversion of polaritons was also dis- 
cussed. The numerical calculations have been performed 
by means of stochastic evolution of a reduced state vector 
in the single particle Hilbcrt space, separately from the 
stochastic environment evolution. The resulting stochastic 
differential equations allows reasonable reduction of the 
problem which by direct numerical integration of the den- 
sity matrix would have being cumbersome. For this matter, 
the stochastic technique turns into a very practical tool for 
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